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Abstract 

The phase diagram for a general model for Cuprates is derived in a mean-field approximation. A 
phase violating time-reversal without breaking translational symmetry is possible when both the 
ionic interactions and the local repulsions are large compared to the energy difference between the 
Cu and O single-particle levels. It ends at a quantum critical point as the hole or electron doping 
is increased. Such a phase is necessarily accompanied by singular forward scattering such that, in 
the stable phase, the density of states at the chemical potential, projected to a particular point 
group symmetry of the lattice is zero producing thereby an anisotropic gap in the single-particle 
spectrum. It is suggested that this phase occupies the "pseudogap" region of the phase diagram of 
the cuprates. The temperature dependence of the single-particle spectra, the density of states, the 
specific heat and the magnetic susceptibility are calculated with rather remarkable correspondence 
with the experimental results. The importance of further direct experimental verification of such 
a phase in resolving the principal issues in the theory of the Cuprate phenomena is pointed out. 
To this end, some predictions are provided. 
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I. INTRODUCTION 



It became quickly apparent following the initial discovery^ that the high superconduct- 
ing transition temperatures in the Cuprates was not the only interesting feature of these 
compounds. The properties in the normal state(s) are unlike other metals and new concepts 
are required to understand them. Moreover the superconductivity can be understood only 
if they are understood since the fluctuations responsible for the superconducting instability 
are the property of the normal state. 

The extensive investigation of the Cuprates 2 -," 3 - has led to a consistent set of experimental 
results from which a universal phase diagram, fig. (1), for these compounds was drawn 4 ^, 
and for which considerable further evidence has been adduce d 3 ! 7 ' 8 ! 9 ' 10 This diagram was 
drawn on the basis of the measured properties, thermodynamic as well as transport, which 
show characterestic changes across the lines drawn in all the cuprates. A fundamental aspect 
of the phase diagram is the existence of a putative quantum Critical point (QCP) inside the 
regime of compositions for superconductivity. Besides the superconducting region, there lie, 
as marked in fig. (1), three distinct regions emanating from the QCP. This paper is concerned 
primarily with region (II) of the phase diagram, but as explained below, its nature is tied 
to the nature of the QCP and the fluctuations above it and understanding it amounts to 
understanding the essentials of the physics of the Cuprate phenomena. 

The QCP was anticipated by the marginal fermi-liquid hypothesis- 1 , which proposed a 
scale invariant spectrum of fluctuations to explain the normal state properties in region 
(I) and led to a set of predictions for the single-particle spectra and transport properties 
which have been verifie d 1214 . The suggestion— that the same spectrum provides a glue for 
superconductive pairing is also consistent with that deduced from approximate inversion of 
angle-resolved photoemission (ARPES) data^ in the superconductive state. The proposed 
spectrum has the form 

7mx(q,w,T) oc u/T, foru^T, (1) 
oc constant, for T <C u> <C u> c , 

where u c is a cut-off, determinable from experiments. The corresponding real part of the 
fluctuations is i?ex(q, uj,T) oc ln(u> c /u>) for u/T ^> 1 and oc ln(u c /T) for u/T <C 1. 

It is characterestic of QCP's in itinerant fermion o 16 ! 17 that a crossover to a fermi-liquid 
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FIG. 1: A schematic of the universal Phase diagram for the Cuprates; the lines shown and the 
crossover occur in all the Cuprate compounds. Other features, small regions of antiferromagnetism 
away from half-filling, charge density modulations, stripes, etc. which occur in one or the other 
Cuprates as well as the spin-glass region are not shown. 

is expected characterized by a line emanating from x = x c . Below another line emanating 
from x = x c , a broken symmetry is expected. These expectations are consistent with the 
experimental phase diagram of fig. (1): characterestic Fermi-liquid properties are found deep 
in region (III) and in region (II) an anisotropic gap (pseudogap) in the one-particle spectra 
at the chemical potential is observed^. Transport- and thermodynamic propertie o 19l2 ° also 
show characteristic changes below the line T p (x). However the nature of region (II) remains in 
doubt. There is plenty of empirical evidence to conclude that region (II) does not universally 
break translational and/or spin- rotational symmetries. If it is a phase of broken symmetry, 
it is of an unusual nature. 

A knowledge of the ground state and excitations of region (II) ties together every essential 
aspect of the Cuprate problem. Given a QCP, region (III) follows automatically. Given the 
success of the predictions following from the phenomeno logical spectrum of Eq.(l), one can 
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be reasonably certain about the functional form of the spectrum. Since the broken symmetry 
in Region (II) is due to the condensation of the fluctuations in region (I), a knowledge of the 
former specifies the nature of the latter. It answers the question, the fluctuations of which 
operator have the spectra given by Eq. (1)? There is experimental evidence^ that Eq.(l) 
also gives the glue for the pairing instability for superconductivity. Therefore a knowledge 
of region (II) specifies the pairing mechanism as well. 

A hint is available about the possible broken symmetry in Region (II) from the Raman 
spectra, which measures the q — > part of the particle-hole fluctuations in even parity point- 
group symmetries. In region (I), it is firmly established that the spectrum is consistent 2 ^* 
with the predicted form, Eq. (1). For conserved quantities like the density or spin-density, 
the correlation functions must be proportional to q 2 at small q. Therefore one infers that 
the observed singularity is in the correlation function of an un-conserved operator in some 
irreducible representation (IR) of the lattice at q — 0. (This IR is inferred from Raman 
experiments to be the B lg representation.) To be consistent with the existence of a QCP, 
the lowered symmetry in region (II) must then be due to the condensation of such an 
operator. In the ordered state, the singularity in the Raman response must then disappear; 
this is also consistent with experiments in region 

An operator satisfying the above requirements is the current operator at q = in the B\ g 
representation of the lattic e) 23 , 24 . In a microscopic theory of the cuprates 4 ^ based on a model 
suggested earlier—, it was predicted that the region (II) violates time-reversal symmetry by 
a ground state which carries currents in the O-Cu-0 plaquettes in each cell. An experiment 
using ARPES with circularly polarized photons was suggested 27 to identify the T-breaking 
phase phase. Some detailed experiments have borne out these predictions 2 - New classes of 
experiments have been suggested to further verify the predicted phas o 29 ' 30 . 

The violation of time-reversal symmetry is shown here to be accompanied by singular 
attractive forward scattering due to the generation of a retarded 1/g 2 interaction. Due to 
such a scattering, the density of states of one-particle spectra at the chemical potential, 
projected to a particular (IR) of the lattice tends to zero. The resulting stable state is 
shown to have an anisotropic gap in the single-particle spectra at the chemical potential^. 

It is worthwhile summarizing the principal results observed in region II of the phase dia- 
gram. The principal property in this region is an anisotropic loss of low energy excitations. 
First discovered in measurements of Knight shifts, this was soon found in the specific heat 
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measurements^, as well as in direct measurements of the single-particle spectra through 
ARPES 18 . Transport properties, Raman scattering and neutron scattering also show char- 
acterestic changes in their temperature (frequency) dependences below T p (x). This region 
is often said to have a pseudogap. This is a misnomer, (which I am forced to persist in 
using.) ARPES measurements show (see for example, figs. (26) and (28) of Ref. (3lll) 
that the anisotropic gap is as robust as the superconducting gap at a corresponding reduced 
temperature. 

The principal purpose of this paper is to give details of the results presented earlier^ and 
to address several related new questions. I focus here on the broken symmetry as well as 
the generated anisotropic gap and the equilibrium properties and the single-particle spectra 
in region (II). The density of states, the specific heat and the magnetic susceptibility are 
calculated with remarkable correspondence with the experiments. The transport properties 
will be studied in the near future. The paper is organized as follows: The motivations for 
the specific model introduced for the Cuprates are reiterated in Sec. II, where the physics 
of time-reversal violation or local circulating currents is related to the screening of charge 
fluctuations special to the cuprates. Time-reversal violating (TRV) states are derived in 
mean-field calculations in Sec. III. A phase diagram similar to fig. (1) is derived. The 
current patterns and the symmetries in this state are also discussed. In Sec. IV, the coupling 
functions of the collective fluctuations towards the TRV states to the fermions are derived 
to show why time-reversal violation must necessarily engender a fermi-surface instability. 
The stable state is shown to have an anisotropic gap at the chemical potential. The single- 
particle spectrum and the thermodynamics are derived and compared with experiments and 
some predictions for further experiments may be found in Sec. VI. 

An important new result is that the derived coupling of the fluctuations in Region (I) 
to the fermions suggests that the phenomenological spectrum of Eq.(l) should be modified 
in one respect: xisi^iT) should be replaced by x(q>k, u,T) where at small momentum 
transfer q, the k dependence near the fermi-surface varies with the dominant IR of the 
fluctuations. 

Various sundry issues, including the limitations of the theory are discussed in the con- 
cluding section. 
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II. THE MINIMAL MODEL FOR CUPRATES 



A microscopic model for the Cuprates must take into account that the phenomena ob- 
served in them hardly ever, if at all, occurs in any other known compound. The effective 
Hamiltonian for the Cuprates must therefore reflect some unique features of their solid- 
state chemistry. Consideration of this issue led to the 3-orbital model with longer-range 
interactions as well as the local repulsions on Cu and on O 25 . 

Cuprates have the unique feature that the ionization energy of the Cu ++ state, is nearly 
equal to the affinity energy of the O state, i.e. E(Cu + ) — E(Cu ++ ) is nearly equal 
to E(0~~) — E(0 ) 26 . The corresponding difference is much larger than the one-particle 
transfer integrals in other transition metal oxides. Therefore the charge fluctuations in the 
metallic state occur almost equally on Copper and on Oxygen ions. The 4s-p orbitals on Cu 
are too far away in energy to screen the mixed- valence fluctuations on Cu intra-atomically; 
the 3s state of O are similarly too far away for intra-atomic screening. (These orbitals, 
when they control the screening, can be eliminated in favor of a renormalized effective 
repulsion parameter, the U of the Hubbard-model.) Charge fluctuations can then only be 
screened by complementary charge fluctuations on the neighboring ions. A consequence is 
that the screening length is close to the nearest neighbor Cu-0 distance giving large effective 
ionic interactions which summed over the nearest neighbors are of the same order as the 
local repulsion energies (U) on the Cu-or O ions. Since both the ionization-affinity energy 
difference and the transfer integrals are smaller than the ionic interactions, there is no small 
parameter in the physics of the Cuprates by which the ionic interactions may be absorbed 
to start with a simpler model such as the Hubbard and t — J model proposed^ 2 - for the 
Cuprates. 

Whether the ionic interactions introduce qualitatively new properties can only be decided 
by calculations using them. I show in this paper that this is possible using approximate but 
systematic methods. For a range of ionic interactions and doping, screening can be achieved 
coherently through occupied wave-functions which are complex admixtures of orthogonal 
orbitals in a unit-cell, i.e. by currents in the elementary O-Cu-0 plaquettes. Such coherent 
screening represents time-reversal violation. 

The minimal model Hamiltonian 25 is on the basis of three orbitals per unit cell, d,p x ,p y : 

H = K + H-n t + Hj® t . (2) 
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K is the kinetic energy operator, 

K = Y^ tdUdka + 2 t pd d^ a (s x (k)p xka + s y {k)pyk a ) - At pp s x (k)s y (k)p^ k{r Py ka + H.C. (3) 

k,cr 

Here a particular choice of the relative phases of the x and y orbitals in the unit cell has been 
made, s X:V (k) = sin(k x a/2, k y a/2) and for later, c X;y (k) = cos(k x a/2, k y a/2) and s 2 xy {k) = 
(sin 2 (k x a/2) + sin 2 (k y a/2)). The local interaction on the Cu and the O orbitals are 

Hint ^ y Ud f^dia ^di~a Up {rip X icr W>p X i—a ~\~ Tlpyj^Tl pyi—a) (4) 

and the nearest neighbor interaction between the Cu and the O orbitals are 

Hint = X! V n di{ni+nn,p x + ^i+nn,p y ), (5) 
i,nn 

where (i + nn) stands for nearest neighbors of Cu in a cell i with (p x ,p y ) orbitals in the 
(x, y) direction. In momentum space, this is 

H-nl =2^ E c x {q) e£ +q(J d ka p^^p^a' + x^y (6) 

kk'q, uu' 

We may also include the nearest neighbor interactions between the orbitals or longer 
range interactions. Such interactions do not change the essential results derived here pro- 
vided they are small compared to V. In the calculations in this paper (renormalised) energy 
difference between the Cu and the O orbitals is taken zero. It is important that in CuO, 
Q ~ 0{t pd ). Taking it as zero simplifies the calculations presented; the principal effect of a 
finite €d will be mentioned. 

In the limit that (Ud,U p ) ^> {t p d,t pp ), a good approximation 3 ^ for small deviations |x| 
from half-filling consists in replacing 3 ^ 

tpd ► tpd ^-pdl^l) tpp * tpp ^pp|^|; (7) 

where \x\ os the deviation from half-filling; x > for holes and x < for electrons. The 
Lagrange parameter appropriate to this transformation is absorbed in the chemical potential. 
A more general and messy mean-field calculation, would consider separately the average 
occupations in the oxygen and copper orbitals and renormalize t p d, t pp accordingly. It is not 
expected to change any of the essential results. 

An alternate treatment 4 to handle the large local and ionic repulsions is to start from the 
strong-coupling limit by diagonalizing a unit-cell exactl y 34 ! 35 by cutting of its kinetic energy 



connection to its neighbors and then to introduce the kinetic energy connecting the lowest 
energy many body states in each cell. This is the generalization of the procedure by which 
the t — J hamiltonian is derived from the Hubbard model. Starting from the model of Eq. 
(1), this leads minimally to an effective SU(2) x SU(2) hamiltonian plus the constrained 
kinetic energy; one £77(2) represents the usual spin-space while the other represents the two 
lowest many-body states obtained in each cell by the first step in the procedure. 

H int = c7(-l/2 + (Ji ■ <Tj)(A + nBrj). (8) 

Here r's are in the orbital space and A, B give the anisotropy of interactions in that space. 
This hamiltonian leads to solutions with the same symmetries as in the approach below. 

III. DERIVATION OF TIME-REVERSAL VIOLATING STATES 

Time-reversal violating (TRV) states are shown here in a mean-field approximation to be 
possible stable states for the model in a range parameters. To construct a mean-field theory, 
write the interactions of Eq.(5) in a separable form. 

#S = -v e 4£t E 4^- (9) 

i=1..4,k,q,(T k',cr' 

A ££l = ( s *( k )ri,k+q j( A >( T ± s 2/ (k)pi k+q (T d ki(T ) (10) 

A< k£c = ( C *( k )rf,k+q j( A >( T ± fy(k)pl,k+q, < A<T) ■ (11) 

The superscripts (1,3) and (2,4) refer to the plus and minus signs respectively. 

A model such as that in Eq. (2) is a severe simplification of the problem. It may be the 
minimum essential model but the parameters that faithfully represent the actual problem 
through it are knowable to factors of 2 at best. Staggered flux phase o^ 1 ^ , charge-density 
waves^, spin-density waves, charge-transfer instabilitie o 39 ' 40 etc. may be produced in mean- 
field by the same Hamiltonian by varying different parameters. It is not the purpose of this 
paper to find the phase diagram for all these instabilities in the space of the parameters of the 
Hamiltonian, but to show that the trial states chosen for investigation by discarding those 
that can be excluded on the basis of experiments occur for reasonable physical parameters, 
and let experiments decide whether they indeed occur. 

Accordingly, I choose to look only for mean-field solutions which do not break transla- 
tional symmetry and/or spin-rotational symmetry. For such solutions, one needs to consider 
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(2) 

only the zero-momentum transfer, q = 0, part of the Hamiltonian H^j.. We can therefore 
introduce possible mean-field order parameters Riexp{i4>i) and rewrite 

#3 = -^£[£4L - 5 q ,o^/Vexp(-i0 i )]E Ag q)ff - ^iyVexp^)] 

q,« k,cr k',(T 

+ £ - i2, £L4» >ff exp(^) + AgoV exp(-i0 i )], 

i k,cr 



^ e x P (^) = ^E(4,U- (12) 

k,<r 



The sum over k in i?'s is restricted to the occupied part and the expectation value is in the 
state to be determined self-consistently. 

The mean-field Hamiltonian with order parameter at q = can only be a 3 x 3 matrix 
in the space of the three orbitals per unit cell at each k. Furthermore, the effective single 
particle energy for the p x and the p y orbitals must be kept degenerate; the renormalization 
of their difference from the energy of the d-orbital can produce a charge transfer instability 40 
which is generically of first order and therefore not of interest. ( We will take the point of 
view that such renormalizations have already been taken into account in the bare parameters 
of the theory.) This leaves only three functions of k in the off-diagonals to be determined. 
The q = part of transforms simply as the existing kinetic energy and cannot lead to a 

(2) 

broken symmetry. If a real order parameter is constructed from the rest of H^J., a structural 
distortion of the unit cell is a necessary accompaniment. We discard this possibility as well. 
So the order parameter must be non-real. As will be obvious below, and A^ lead to 

(2,) 

the different domains of the same phase. So a mean- field theory of H- n l of Eq. (5), not in 
conflict with experiments, can produce only two distinct order parameters: 

RjeM^i) = ^E(41.> = ±^E[^fe + k,A.) - *y<pk,A>]- (13) 

k,<r k,cr 

RnexpWn) = VJ2(A^) = ±V Y\ C M*,J^) ± ^feA,.)]- (14) 



In the mean-field approximation the expectation value of the constrained kinetic energy 

(2) 

plus H> nt is minimized to determine the order parameter R. For the state 0/, the mean-field 
hamiltonian is 



H I mf = K + R]/V - R 1 expH0 7 ) Y>*toPX4** ~ *vWP&«A*A + H C - (I 5 ) 

k,ir 
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For the state 0//, the mean-field Hamiltonian is 



H% f = K + R 2 H /V - R n expi-i^j) Y\c x {k)p^ i(T d^ a - c y (k)p^ a d k , a ] + H.C. (16) 

k,cr 

For both states the first and the third terms of the mean-field Hamiltonian written in the 
space of the orbitals dk,p x k,p y k is the matrix 

/ 

Hmf = 



2i pd e ie *Ms x (k) 2i pd e ie y^s y (k) ^ 

2i pd e- i9 *Ms x {k) 4t ppS ,(k)s y (k) 



\ 2i pd e- i0 yM Sy {k) Ai pp s x {k)s y {k) 



For the state 9/ 



tan 6^ 
tan 6. 



Ri sin . 



2t pd + Ri cos <pi ' 

— Rj sin <pj 
2t pd — Ri cos (pi 



while for the state Qn 



tan(e x (k)) -- 
tau(0 w (A;)) = : 
The unitary (gauge) transformation 



R n c x sin (fin 



2t pd s x + R n c x cos 4>n 
RnCy sin (j)n 
~2i pd s y ± RnCy cos <pi i ' 



dk,a -> d K(7 ; p xXa -> p X)kjCr exp(i# x ) ; p yikjCT -> p W|ki(T exp(i0 s 



on Eq. (17), transforms it to 

( 2i pd s x (k) 
2i pd s x (k) 



(17) 



(18) 



(19) 



2t pd Sy (k) 

4^e'(«--«-) ai (k) ay (k) • (20) 
{2i pdS y{k) Ai pp e^y~ e ^s x {k)s y ik) , 

Note that 

0(k) = (0«(k)-0 y (k)), (21) 

is the invariant /Zua; associated with the elementary O-Cu-0 plaquette. (The correspondence 
of the theory here to the generation of time-reversal violation in SU(3) model of elementary 
particles through the Cabibo, Kobayashi, Masakawa SU(3) mass matrix^ may be noted. 
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The time-reversal violation there is in an abstract space; here the phenomena is more vividly 
pictured.) The problem therefore reduces to determining this invariant by minimization of 
the Free-energy. 

Returning to the minimzation, since at any k the trace over the three bands is conserved, 
the change in ground state energy can be calculated from the change in energy of the unfilled 
part of the conduction band alone. In an Appendix these are calculated to leading order in 
tpp/tpd- The value of the phase <piji at the minimum is found to be at 

= ±tt/2. (22) 

With e<2 = 0, and the deviation from half-filling \x\ <C 1, 9i$ ^ at T = 0, only if 

V Z 2a 1)2 -^^/4, (23) 

where are numerical coefficient of order unity, which depend on the details of the band- 
structure and whether the state is 0/ or Qjj. The equality in Eq.(23) provides the hole 
(x c > 0) and the electron (x c < 0) densities at the QCP's. In Appendix A it is shown that 
the simplest model favors the 0// phase. Only the possibility of the 0/ phase was considered 
earlier 5 . (23) gives that for V ~ t p d — 2t pp , x c ~ 0.1 which, in view of the approximations 
made, is surely a stroke of luck. 

In the Appendix the transition temperature is also derived to find 

T g oc Ef/ln(l — x/x c ), x < x c , (24) 

where Ef is the Fermi-energy in the conduction band. At x — > 0, this result is an artifact 
of the model with taken zero. (At x — > 0, the effective masses in all bands are infinite, 
so that for = the bands are degenerate and the transition occurs for V — > at T — ► 
or infinite T at V finite.) Moreover competition of antiferomagnetism and the TRV state 
has not been considered here. The result is therefore only valid for x close to x c <C 1. Eqs. 
(23,24) give the correct general shape for the phase (II) marked in fig. (1) for x below near 
x c . 

The mean- field wave-function is made up of products of |k, up to the Fermi- vector. 
In the original gauge, i.e. in which Eqs. (2) and (15) are written, the new (TRV) conduction 
band wave-functions for either state 0/ or 0/ to leading order in t pp /t p d are, ( the spin-index 
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is suppressed), 



<. 1 |0) = (iV fc )- 1 [a+ + 2 




((si - s 2 y ) cos(0 (1 , 2) (k)) + ^4sin(^ (1)2) (k))<)]|0) 




(25) 



e£ = 2t pd s X j / (k) and TV^ are normalization factors, 0i(k) are variationally determined: 



9 2 (k) = ±(< Rn > /2t pd )(cot(k x a) ± cot(k J/ a)), for state Q n . 

These expressions are obtained from (18) and (19) for < R > /t pd « 1 and k X) k y not close 
to zero for 6//. 

A. Symmetries of the Mean-field TRV Hamiltonian 

It is useful to recall that a complex hamiltonian which cannot be transformed to a real 
operator by any unitary (gauge) transformation does not commute with the time-reversal 
operator 7Z and therefore has eigenstates which break time-reversal invariance. Correspond- 
ingly, its eigenstates cannot be transformed real by any unitary transformation. This is 
true of the wave-functions of Eq. (25,26). So Time-reversal is violated in the phases found. 
Further, if 



72. | k, 9, a) is not an eigenstate of H(9) but is of 7Z^ 1 H(9)7Z with the same eigenvalue e(k, 9). 

Consider the state Oj. Physically, 9i/2 and —9i/2 are the phase difference between the 
d-orbital and the p x and p y orbitals respectively in each unit cell. The total phase difference 
9\ in going around any of the four triangular plaquettes in each unit cell is gauge invariant 
and can be distributed between the three legs of the plaquettes in any way one chooses by 
an appropriate gauge transformation. If we include additional interactions which produce 
additional phase shifts between the vertices of the plaquettes, (for example a repulsion 
between electrons on the neighboring oxygen orbitals produces in a mean-field approximation 
a phase shift between the p x and the p y orbitals of each unit cell) it can be absorbed in the 



0i (k) = ± < Ri > /2t pd , for state 6/, 



(26) 




(27) 
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invariant flux. A finite flux in a plaquette is equivalent to a current circulating around the 
plaquette. The current pattern for (one of the two domains of ) the state 0/ is shown in 
fig. (2). This current pattern can be found by transforming the wave-functions to real space 
and calculating the expectation value of the current operator between the sites of the lattice. 
The magnitude of the current is of 0{6it pp t 2 pd a' i ). 

The symmetry of the current pattern follows from the properties upon applying the 
reflection operator about a mirror plane rh on the wavefunctions (25,26). When rh is the 
x = or the y = plane, the wavefunctions are eigenstates of cr^- But reflection symmetry 
is broken about the crystalline mirror-planes x = ±y. Note that if we consider properties 
which depend on the modulus of the wavefunctions, for example the charge density, mirror 
plane symmetries existing for 8 = are preserved for 9 ^ 0. A current sensitive experiment, 
(sensitive linearly to the wavefunction and its gradient) can therefore detect a Time-reversal 
breaking phase through the difference in reflection symmetries but experiments sensitive 
only to charge-distributions cannot. 

It may also be seen from the wave-functions that four-fold rotation (C4) about an axis 
through the Cu's is not a symmetry but the rotation followed by time-reversal C4JI is a 
symmetry. Also the center of inversion is preserved. The simplest lattice symmetry for 
the Cu-0 compounds is 4/mmm; the lowered symmetry generated below the transition is 
4/mmm in the double-group notation^- 



FIG. 2: (a): The current pattern in the time reversal violating state 0/ and (b); the current pattern 
in the time reversal violating state 

In the state 0//, reflection symmetry is lost about the crystalline mirror-planes x = 
and y = 0. The states with plus signs between the p x and the p y parts in the expression 
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for 62 keep reflection symmetry about x = y mirror planes but not about the x = —y 
mirror planes; the converse is true for the state with the minus sign. The current pattern 
for the former is also shown in fig. ( 2). Inversion is lost in this phase but the product of 
inversion and time-reversal is preserved. Such phases are called magnetoelectrie^. Starting 
from lattice of symmetry 4/mmm, the reduced new symmetry is mmm. These issues are 
further discussed in Ref.(— ), where extensions to more complicated lattice structures is also 
given. Note that the 0// phase has four domains. 

One may check^ 3 - that neither the 0/, i.e. 4/mmm nor the 0// phase, i.e. the mmm 
is piezomagnetic; the application of a uniform magnetic field does not linearly change the 
symmetry of the unit-cell. 

The polarized ARPES experiments^ are consistent with the state 0// in region II of the 
phase diagram in fig. (I) 44 . It follows that if the superconducting transition is a continuous 
transition, the symmetry breaking of the pseudogap state, both time-reversal and inversion, 
should continue in the superconducting regio n 30 ' 45 to the left of the QCP in fig. (1). An 
experimental evidence for this would constitute a further direct test of the theory of the 
Cuprates presented here. 

The interesting property that chiral spin correlations accompany orbital currents of 
fermions ought to be mentioned 46 . For three sites (i,j,k), and current operators Jjj from 
site i to site j, etc., the expectation value 

(JijJjk-Jki) oc (Sj • (a, x B fc )), (28) 

where Sj is the spin-operator at the site i. (The order parameter 9 is proportional to this 
current expectation value with (i,j, k) forming the nearest neighbor O — Cu — O triangular 
plaquettes.) Therefore chiral spin order always accompanies the orbital current order. 

IV. THE FERMI-SURFACE INSTABILITY OF THE TRV STATES 

The TRV instability derived above is due to the mixing of states from different bands 
(equivalently due to mixing of different correlated orbitals in a unit-cell if we start from the 
strong-coupling end) due to large enough screening interactions. It is not a Fermi-surface 
instability and by itself does not create a gap at the chemical potential. I show here that such 
a TRV state cannot exist with a normal Fermi-surface. The stable state has an anisotropic 
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gap at the chemical potential. TRV as well as translational symmetry is preserved in the 
new state. 

Consider the classical fluctuation regime in the transition to the TRV breaking phase. 
(This should be distinguished from the regime I of the phase diagram in fig. (1)) The 
simplest representation of the current fluctuations in this regime has a propagator of the 
form: 

X c (q,cu) oc {iu/ lQ + K 2 \q\ 2 + e)~ x (29) 

e is the distance to the transition \T g — T\/T g but the region close to the QCP at T g = is 
excluded from the considerations here. Since the fluctuations are of an unconserved quantity, 
the leading term in their damping is a constant denoted here by 70. The dynamics belongs 
to Model A in the classification of Halperin and Hohenberg 47 . Since the order parameter 
is also a discrete variable, the fluctuations have a finite coupling to the fermions in the 
long wavelength limit. This is shown explicitly below. Due to the infinite range of the 
fluctuations, the normal Fermi-surface cannot be stable. This is seen in a calculation of 
the Fermion self-energy due to exchange of the fluctuations of Eq. (29). One finds that for 
|k- k f | -> 0, 

E(M) = |7(M| 2 E r du ImXc ^^ . r (coth(o;/2r)+tanh(e k+q /2r))- (30) 

q J -00 V + U) — 6 k+q + 10 

g(kf) is the coupling function of the fermions to the long wavelength fluctuations, derived 
below. For a classical transition, we may take u <C T and consider the pole in \ c to be 
much more important than the denominator in Eq. (30). Then on removing u from the 
denominator and noting that the integral of lm\ c over u is 0, one is left to evaluate only 

S(k, ,) » 2T h{ k f)l r S - _ J q - - s £ ^ - ^ - f)2 ■ (31) 

This gives for small v and small k — kf that 

£(k,i/) oc \l(kf)\) 2 {v + |k — kf | + e)" 1 hx{v + |k — kf | + e) in 2 dimensions, (32) 
oc \^{kf)\) 2 h\{y + |k — k f | + e) in 3 dimensions. (33) 

Recall that the self-energy due to (unscreened) Coulomb interactions in the Hartree-Fock 
approximation has a In |k — kf | divergence in 3d and a |k — kf | _1 divergence in 2d. Here 
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because of the retarded nature of the effective "Coulomb" interaction, the singularities in 
|k — kf| and v are identical. The Hartree-Fock singularities are eliminated when screened 
interactions are considered. In the present problem, such singularities are protected in the 
vicinity of the TRV transition. 

The renormalized fermi-velocity is given by 

v(k / )/v /0 = z{\ + vf Q ■ 0E/0q)| k/lMI (34) 

where vy is the bare fermi-velocity and z is the quasi-particle renormalization, z — (1 + 
dH/du) -1 ^^^. Therefore the fermi-velocity is not renormalized due to the cancellation in 
the divergence in the frequency and momentum dependence of the self-energy. However the 
quasiparticle density of states at the fermi-surface is given by 

N(k f )/N = (1 + vjo 1 • dZ/dq)- 1 ^,,. (35) 

Therefore the partial density of states at the fermi-surface in the irreducible representation 
(IR) of |7(&/)|) 2 tends to as e — > 0. If the partial density of states in any IR approaches 0, a 
compressibility in that IR approaches signifying an instability. The cure to this instability 
is shown through an explicit calculation below to be a state with a gap in the density 
of states in the IR in which g(kj) is maximum. The connection of this situation to the 
Landau- Pomeranchuk instabilities is discussed in Ref . ()49l ) . Since the instability necessarily 
accompanies the classical fluctuation regime of the TRV transition, the two must occur 
together, as is also shown below. 



A. Coupling of the Fluctuations to the Fermions 

In this section the coupling functions 7(k) are derived in the fluctuation regime to the 
TRV state and the TRV state itself. To this end, the mean-field hamiltonian, Eq. (15,16) 
is supplemented by the terms giving the fluctuations in Eq. (12). 

In the 0/ phase, the fluctuations either change the magnitude of the order parameter: 
[9i — > 9\ + 5# q ) or admix the states of the other domain (<p — ► ±7r/2 + <50 q ). The last two 
terms of Eq. (12) have been used in H m f. The first term is used to generate the fluctuations. 
To this end, one writes 

v4k, q = (#i + d&q/V) exp[i(0 o + ^0q)] + Incoherent part. (36) 

k,q 
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The incoherent part is the fermionic part of the fluctuations which is expressed in terms of 
the new particle-hole states in the 6 phases. The coherent part are the bosonic fluctuations. 
In the fluctuation regime to the 6 states, 9 = and the only fluctuations are the amplitude 
fluctuations <5# q in the simple small fluctuation approximation used here. 

The Fermion-Boson coupling hamiltonian follows directly by finding the matrix elements 
of the fermionic part of A kq in the wavefunctions given in Eqs. (25,26). The results are 

Hfb = £ cl,eCK0 (7p(k, k')[<%_ k0 + 50L (k _ kO ] + 7a(k, k')[<%_ k0 + ^l (k _ k/) ]) (37) 
k , k ' 

In the next section we will be interested in the coupling only in the forward scattering limit 
(k — k' — > 0). In this limit, for the Oi phase, 

7a(k, k) fa VQ\ (1 + cos(k x a) + cos(k y a)) , 7 P (k, k) fa V9\ (cos(k x a) — cos(k y a)) . (38) 

The energy of the phase mode at long wavelengths is calculated in Appendix A to be 

~ %l (39) 

The energy of the amplitude mode can be shown to be higher by a numerical factor of 0(2). 

In the fluctuation regime, the only coupling is to the amplitude mode. The leading 
forward scattering coupling has the dependence, 

7/ toct (k, k) ~ V/2 (cos(k x a) - cos(k y a)) . (40) 

In the O// phase the fluctuations may change the magnitude of the order parameter 
(02 — > 2 + 5# q , admix the states of the other time-reversed phase — > ±7r/2 + 50 q ,or admix 
the states of the other spatial domain, i.e. admix states with reflection symmetry about the 
(x + y) if the mean-field state chosen is the one with reflection symmetry about (x — y) 
or vice-versa. We will call the operator for the last 5^ m>q . So in addition to the two kinds 
of coupling in Eq. (37), a coupling to fluctuations S9 mj(l + 59^^ with a coupling function 
7 m (k + q, k) must be included. 

To ensure a gauge-invariant calculation in the 6// phase, a hamiltonian with a larger 
symmetry than that used to generate the mean-field state must be used. A minimum such 
Hamiltonian includes both the A\A 3 and A\A^ terms of Eq. (9). We can then follow the 
familiar procedure to generate the coupling functions in the forward scattering limit for 
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the three modes already described: For the amplitude mode, phase mode and the "mixing 
mode" we get, respectively 

7o(k, k) ps V9 2 (cos(k x a) + cos(k y a)) (41) 
7 P (k, k) ps V9 2 (sm(k x a) ± sm(k y a)) (42) 
7 m (k, k) Ps V^cos^a) — cos(&; y a)). (43) 

The energy of these modes in the long wavelength limit are estimated in the Appendix. 
The energy of the phase mode Q<p is about the same as that in the phase 0/. the energy of 
the mixing mode at long wavelengths Q m is estimated in the appendix to be about 1/2 of 
these. Therefore the important mode to consider for stability in the state 0// is the mixing 
mode. 

The forward scattering with the amplitude fluctuations in the fluctuation regime has a 
dependence, 

7/; uc t(k, k) ps V/2(s\n(k x a) ± sin^a)) (44) 

I have only considered small fluctuations in the above. The model also has interesting 
topological fluctuations (mentioned in the concluding section.) 

B. Anisotropic Gap (Pseudogap) in the TRV state 

In the first part of this section, I presented an argument from T > T g that the TRV 
phase has a tendency to develop an anisotropic gap at T = T g . Now, I present a calculation 
coming to the same conclusion from T < T g . In effect, the zero order collective modes in 
the TRV phase, whose frequencies are estimated above, are shown to be unstable due to 
renormalizations by Hfb with a normal fermi-surface. 

The effective Hamiltonian in the TRV phase is 

H t?f = E e ( k - *)cL<V + E ^(q)6 + (q)%) + H FB . (45) 
ko- q 

Hfb is given by Eq.(37). The annihilation, creation operators of the fluctuation modes are 
denoted by b,b + . The fluctuations with the lowest energies are of interest. They are the 
fluctuations of the phase for the 0/ phase and the "mixing" -fluctuations for the 0// phase. 
The operators 6(q), b + (q) will only refer to these. We may eliminate these fluctuations and 
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Hfb to generate a retarded effective interaction between the Fermions of the form 

Yl 7( k > k + q)7(k' - q, k')x c (?, w)4 +q(J c kCT c^_ qCT ,c k . CT /. (46) 

k,k',q,cr,(j' 

where x c (g, cj) is the propagator for the relevant fluctuations. Xc m &y be approximated by 
— with a cut-off such that e(k'),e(k' + g) are both within about u c of the chemical 
potential. The effective interaction is repulsive for states at higher energies. We use only 
the forward scattering part of the interactions, i.e. in the limit of small q,and for k, k' close 
to the Fermi-surface: 

g(k, k') = -Lim^ 2 7 (k, k + q) 7 (k' - q, k')/fi q , (47) 

with a cut-off of order u c ~ f2 . If Sl — ^ 0, the cut-off u c is approximately the larger of Qq 
or the damping of the collective mode 70 |e of Eq.(29). 

The instability seen through Eq. (30) may also be seen through the equation of motion for 
the fluctuations < c^ +q)0 .Ck l0 - >. Just as for Landau- Pomeranchuk instabilities, no solution 
for such equations can be found for real frequencies in the limit of small q 48,4 «. The cure 
to this instability is a new state in which particles are annihilated and created by operators 

CT , Ck,(T respectively, such that a stable homogeneous equation for 

)(q,u). (48) 

at v = can be found at q. Correspondingly, in the stable state, the one-particle eigenvalues 
are renormalised from €k6»o- to E^ a . Actually, the interactions 7(k, k') are renormalized also. 
But such strong-coupling corrections are not considered here. The particular form for 
and the relation of c's to c's, etc. is dictated by the self-consistency condition derived below. 

A stable state is chosen in which < c£ + >, and < 6 q >, < 6 q > acquire non-zero 

values for small q. ( In taking the limit, k + q and k should be on opposite sides of the 
Fermi-surface; the limit q — ► of < c£ +q Ck i<7 > is not the particle density at k.) The self- 
consistency equations are derived in a mean-field approximation as follows: By a mean-field 
decomposition and minimization of the effective hamiltonian, Eq. (44), 

' 7(k^k + q) 

(K) = ~ ( C k+qaCka) (49) 

kcr Q 

The restriction on the sum indicates that the sum is over states within uo c of the chemical 
potential and that the states k, k + q refer to (occupied, unoccupied) states or vice-versa. 
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The effective mean-field hamiltonian for the electrons is now 

H eff = e k4 CT C kc7 + V ( k > k + q) C k+qaCka = E E k C^ a C ka (50) 

kcr k,q<7 kcr 

V(k,k + q ) = 7 (k,k + q)((6 q ) + (6l q )). (51) 
One may write i? e // in real space as a tight binding Hamiltonian on the lattice, 

He ff = £ f(Ri,Ri)<&c jVT , (52) 

which shows effective long-range, angle- dependent hopping, which is to be determined self- 
consistently. 

The one-particle Green's function for H e ff in the self-consistent Brillouin-Wigner ap- 
proximation has a self-energy which is an infinite continued fraction which may be summed 
formally to give the one-particle propagator 

G(k,u) = 1 —— -, (53) 

v ; w - e k - £>(k, w) ' v ; 

where 

D M = ± lViKk : ( f , (54) 

On the energy shell, at 

cu = E k = e k + D(k,E k ), (55) 

D(k, E k ) = ± (f^)-f(E k ))\V(k,k + q) ^ (56) 

q -^k+q — E k 

This formal solution to (50) is exact in the limit that the number of states (k + q) coupled 
to a given state k is very large compared to 1. It also follows by taking the expectation 
value of (50) in the new state that 

D(k,E k ) = £\/(k,k + q)(c+ +q c k >. (57) 
q 

Comparing Eqs. (56) and (57) and using Eq.(51), 

/+ v 7 (k,k + q)(/(E k+q ) - /(£ k ))((&q) + (gg) _ 
(Ck+qCk) = p p • (58) 

-^k+q — -^k 
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Inserting (& q ),(&t q ) from Eq. (49) and after a simple manipulation I obtain the self- 
consistency condition : 

V(k + q/2, k - q/2) = £ g(k, ^ y'^f ^ + q/2, k' - q/2), (59) 

In Eqs. (59), the sum is restricted to states k' such that e(k') is within about Qq of the 
chemical potential //, as required by the retarded nature of the Fermion-Boson interaction. 

Note that we have already proven in Eq. (55) that the quasi-particle energies in the new 
state obey, 

E ka = e ka + D ka . (60) 

There is no self-energy in quadrature as in the BCS theory or the theory of CDW/SDW. 
This of-course has pronounced effects on the spectroscopic and thermodynamic properties 
in the new state. 



C. Solution near T g 

The self-consistent solution of Eqs. (56), (58) and (59) is a complicated problem which 
however simplifies in the " Ginzburg-Landau" regime which provides the transition temper- 
ature and the leading temperature dependence of the gap below it. In this regime only very 
small q old states are admixed to a given k in the new states. I will assume that the solution 
obtained continues qualitatively to lower temperatures. 

In this regime, note first that since g(k, k') is proportional to 7(k)7(k') for q — > 0, 
V(k,k + q) ~ 7(k) also. Second, E(k) « e(k) for \e(k) — fx\ ^ Qq. The self-consistency 
condition (59) may therefore be rewritten in the limit q — > as 

Here E^,E^ are the new one particle eigenvalues above and below the chemical potential 
respectively. To determine E k , one must satisfy Eq. (60) in a manner continuous with the 
approach to the instability, where <9£/<9k diverges. Accordingly, I seek a solution with 

E> = e k + D(k f ), 
E< = e k -D(k f ), 
(£>'< - e k ) - far |E>'< - fi\ » Q . (62) 
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with the additional condition that D(kf) > for all kf. 

This manner of opening a gap at the chemical potential appears is new. The physics 
of this gap is distinct from that of the gaps of Bloch- Wilson one-electron theory, the BCS 
gap, the Mott gap or the charge- density (CDW)or spin-density wave (SDW) gaps. In the 
mean-field approximation the gap arises due to an effective infinite range hopping exhibited 
in Eq. (52). 

The integral equation (61) for D(k) may be rewritten as 

l = -—f ^i hMl r dep (A tanhimP)) 

n J 2tt D(k f ) Jo m } l + cosh(Pe)sech(PD(kf)' 

At temperatures for which f3D(kf) <C 1, Eq. (63) reduces to 

1 - tanh((3u c /2) j ^P(^\l(k)\ 2 = (64) 
-(3 2 /l2(tanh{(3u c /2) -tanh 3 (/3u c /2)) ±-p(0) J ^(k^D^kf). 

This constitutes the microscopic derivation of a free-energy of the Landau form for the 
instability. 

It is interesting to make connection with the forward scattering parameters of Landau 
theory. This may be done only when the cut-off u c satisfies [3il » 1. In that case define, 

1 r rikx 

"^ P(0) /^f l7(A;/)|2 = F//(2£+1) ' (65) 

where F/ is the spin-symmetric Landau parameter in the £-th " angular momentum" channel. 
The condition for the instability is that the left side of Eq. (64) be zero. This happens at 
a temperature Te below which stability is achieved with a finite value for the right side, i.e. 
with D(k,T) 7^ 0. The transition temperature in that case is 

T e « , Uc ^. (66) 

In Um^r 



In our problem we must take the limit, f3u c 1 for T just below T g , where Qq — > 0. 
Then the cut-off uj c ~ 7o|e|, defined in Eq.(29). The transition temperature T' g is then given 
by 

22 



Here the relationship, Eq.(39) between the collective mode frequencies f2 and the order 
parameter 9 has been used. Eqn.(67) gives T' g m T g for (p Q V 2: t p( fjo) / {Tgt^p) > 1. For 
7o ~ Po 1 an( i x ~ 0-1) this ratio is of O(10 4 ). One should expect that in a better theory, the 
TRV instability is accompanied exactly by the instability to the anisotropic gapped state. 
As explained later, the mean-field theory of the TRV transition is expected to be strongly 
modified by fluctuations. The results for the region near the transition obtained here should 
only be regarded as indicative. 

The value of the gap for T < T g is also given by Eq.(64): 

The angular dependence of D(k) is to be found from the equation (68), subject to the 
condition that D 2 (k) > for all k. A family of solutions to (68) exist, the one that must 
be chosen is that which minimizes the energy by considering the fourth order terms in the 
free-energy. This is given by D(k) which has the same dependence as the kernel \g(kf)\ 2 , 

Z^)oc| 7 (MI 2 - (69) 

Referring to Eqs. (38,40), in the 6/ phase, \g(kf)\ 2 oc (cos(kf x a) — cos(kf y a)) 2 . (This is also 
true in the fluctuation regime to the However in the fluctuation regime to the 6// phase, 
using Eq. (44), \g(k)\ 2 oc sin 2 (kf x a/2 ± kf y a/2), the two signs refer to the two different 
reflection symmetry violating domains of 6//. 

In the 9// phase, the most unstable mode is the mixing mode for which the coupling 
function again has "d-wave symmetry' and the gap is expected to have the shape (cos(kf x ) — 
cosCkfy)) 2 - 

The value of the gap is expected to grow for lower temperatures approximately as in 
Eq.(68). The amplitude of the gap in the one particle spectrum at the chemical potential 
D at T = is w VQT g {x). 

One should ask at this point how the transition to the TRV symmetry is affected by the 
pseudogap appearing simultaneously with it. Recalling that the transition to that phase 
depends in admixing states between the bonding and the non-bonding bands, the effect 
is only of order D /t p d « 1. It would of-course be desirable to derive the TRV plus the 
anisotropic gap in a single mean-field ansatz rather than with the procedure used here. 
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V. SOME PROPERTIES OF THE TRV-PSEUDOGAP STATE 

In (62), the dispersion is determined for /3D(k) <C 1. To calculate properties at 
intermediate and low temperatures we need for (3D(k) > 1. To meet the requirement 
that E k — > ek for \E^ — ji\ ^> D{k), we may choose 

£ k = e k + D(£;)/(l + (e k /e c ) 2 ). (70) 

The derivation of e c is only possible through a solution of the self-consistency equations 
at low temperatures together with the requirement that the number of states below the 
chemical potential is constant. One expects e c to be O(D ). The dispersion in Eq. (70) is 
represented in fig. (3). 




FIG. 3: The one-particle dispersion measurable by ARPES expected for the pseudogap phase. e c 
has been assumed to be = Dq, the maximum gap. The discontinuity D is a function of angle 
around the Fermi-surface. 
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The density of the state in the pseduogap phase may be calculated as 



N p {u) = Y,(l/ir)6(u-Ek). 



(71) 



k 



If e c were 3> -D , this can be evaluated analytically to give 



N p (u) = 2^(0) — arcsin 



( 



D 



CO 



1/2 



) 





(72) 



This increases as \u/D\ ' for \u/D\ « 1. For a general e c , the asymptotic low energy 
density of states keeps this form. For e c of O(D ), the density of states can only be evaluated 
numerically. 

This density of states ignores both the effect of impurity scattering as well as lifetimes 
due to inelastic scattering which are a function of energy and temperature. The latter are 
included below in calculating the single-particle spectral function. 

A. Single-particle Spectral Function 

1. Inelastic Scattering 

To obtain the spectral function A(k,u) measured by ARPES one needs besides Eq.(70), 
the self-energy of the single particle states. To calculate the scattering rate, we first estimate 
the polarizability Xo(k, q, uj) from a single-particle- hole bubble formed from states k and 
k + q. 



We are interested especially in the low energy contributions to imx (q, ku;). These come 
when the initial state k is near a zero of the pseudogap. With such an initial state and a 
the gap in the single-particle spectrum D(q), the integrals in Eq.(73) may be evaluated to 
yield the result that 



7m X o(k, q, u) oc ]T 5 (uj - (£ k+q - £ k )) [/(^k +q ) - 



(73) 



k 



Imxo(K q, u) ocuQiuj- D(q))/v 2 f 



(74) 
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with a function of the angle q multiplying (74) which has a weak dependence, except when 
k + q is also near one of the zeroes of the gap function. In that case 

Jmxo(k, q, a;) ococ uo / q 2 Q(uo — -D(q)). (75) 

Given that the spectral weight of fluctuations for most (k, q) is below an energy D(q) 
and increases only linearly beyond it, no important modifications are expected due to inter- 
actions. For uj ^> D(q), Imxo(k, q,w) must have the same form as in the normal state for 
all k, i.e. of the marginal Fermi-liquid hypothesis. 

From the estimates above one can calculate the form of the self-energy S(k, uj) for k in 
the vicinity of the region of zero gap. For this calculation the region in which Imx is given 
by (75) is ignorably small. Using (74), a simple calculation gives for uj/D <C 1 that 



IrnL(k,u) oc N(0)u 2 /D . (76) 

For uj/Dq ^> 1, the self-energy must revert to the MFL form. 

For a k away from the zero of the gap, and near energy D(k), no decay through an 
inelastic process is possible kinematically for uj less than D(k). (This is the process where the 
intermediate Fermion state is near the zero of the pseudogap and the scattered polarizability 
makes up the momentum difference.) Therefore for (uj,T) « D(k) the self-energy is 
exponentially small. For (ou,T) > -D(k) it returns to the value Ti n (uj,q,T) without the 
pseudo-gap. So 

IrriE(u, k, T) » sech ( ^2+^)1/2 ) k, T) . (77) 

The spectral function at the wave-vectors away from the zeroes of the gap E(kf) = 
±D(kf) calculated using Eq. (77) and the marginal Fermi-liquid form for S n is plotted in 
fig. (4) for a few temperatures. A 'pseudogap' in the direction kf appears below T m D{kf). 
At zero temperature the Fermi-surface is composed of four points. But at finite temperature 
a gap can be discerned in a direction k/ only if T is less than about D(kf,T). This produces 
the illusion of 'Fermi-arcs', which shrink as the temperature decreases. 
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FIG. 4: Calculated Spectral function in the pseudogap phase excluding linewidth due to impurity 
scattering. The spectral function in any given direction k has been normalized to the value of 
the gap in that direction Dk and the energy is also normalized to the value of the gap. The 
temperatures shown are given by T = nD(kf)/4ir for n=l,2,..5. 

2. Impurity Scattering 

It is important in the Cuprates (and other layered materials) to distinguish the effect 
of impurities in the planes which are expected to acts as point or large-angle scatterers 
from that of impurities in between the planes which scatter electrons in the planes only 
through small angles on the Fermi-surface^ 7 . In the latter case, the scattering rate at a 
point k (calculated in the Born approximation) depends on the local density of states at 
the point, i.e. on the spectral function A(k,u)). This idea is important to understand the 
difference in the scattering rate deduced from single-particle spectra from that in transport, 
the insensitivity of T c to small angle scattering, the Hall effect and the lineshape as a function 
of u and its variation with kf of the spectral function A(k, u) in the superconducting stated. 

Given the form of the spectral function for the pseudogap neglecting impurity scattering 
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in fig. (4), small angle impurity scattering also has interesting effects on the observable spec- 
tral function in the pseudogap state. It follows that elastic scattering will be anisotropic, 
increasing from the nodal region where the spectral function is small to the antinodal re- 
gion, where it is large. It also follows that the small angle scattering rate will be small 
in the antinodal region in the tail of A(h, uS) towards the chemical potential and large on 
the other side. Such features have been observed in experiments^. The detailed compar- 
ison of experimental A(k, u) with calculations should incorporate these effects of impurity 
scattering. 

In the classical fluctuation regime I, above the pseduogap regime, the single-particle 
spectra should also show the frequency dependence calculated in Eqs. (31,32). Since, it is 
only weakly dependent on z/, \k— kf \ and e, it may appear within the experimental uncertainty 
as enhanced elastic scattering, which is angle dependent oc \l(kf)\ 2 . 

VI. COMPARISON WITH MEASURED PROPERTIES OF THE PSEUDOGAP 
PHASE 

The effect of impurities is much less in evidence when one calculates the density of state 
N P (E) by integrating the spectral function over k. The calculated N P (E), for any x at several 
temperatures below T p (x) is shown in fig. (5). To compare, the measured density of states 
by STM is shown in fig. (6). I do not have anything to say about the general approximately 
linear variation of the density of states as a function of energy over the whole energy range, 
which is found both above and below T p . If one looks at the features brought about by the 
pseudogap, good correspondence may be discerned in both the energy and the temperature 
dependence. Of special importance is to note the absence of the BCS type singularity in the 
density of states in the pseudogap both in the calculations and in the experiments in the 
pseudogap regime. In the experiments it appears as one cools below T c (x). (Such singularity 
is present in the pseudogap region in theories based on any translational symmetry breaking, 
for example a charge density wave or a D-density wave.) The variation of the pseudogap 
feature as a function of temperature in the calculated curves also bears close correspondence 
to the experiments. 

The specific heat coefficient, 7(T) = C v /T and the magnetic susceptibility x(T) at differ- 
ent x, calculated from the density of states of fig. (5) are presented in fig. (7) and (8) together 
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FIG. 5: Calculated density of states in the pseudogap phase at several temperatures. Dq is chosen 
2.5T g 

with the experimental results. The specific heat results in Ref. (3) for various x are rescaled 
with a T p (x) to lie as shown in fig. (7) with a common ratio Dq(x)/T p (x) = 2.5. The values 
of T p (x) are within the uncertainty the same as the values where the pseudogap is seen in 
the specific heat or the resistivity. Considering the errors in extracting the electronic specific 
heat from the measured specific heat, no more sophisticated fitting is warranted, although 
clearly the fitting could be further improved by choosing a slightly variying D /T p (x) for 
different x. 

The temperature dependence of the measured x{T) for a particular x is compared with 
experiment in fig. (8) with again D /T p = 2.5. Within the uncertainties in the determination, 
the experimental C V (T, x) /Tx{T, x) is independent of temperature^, in agreement with the 
result here. The agreement of the calculated thermodyanmics with the experiments should 
give one confidence in the underlying theoretical ideas. 
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FIG. 6: Measured density of states at temperatures in the pseudogap phase obtained by Renner 
et a!., Ref. jf^m Bi-2212, by scanning tunneling microscopy 



From Eq.(72), it follows that the low temperature specific heat C v is predicted to be 
~ T 3 / 2 and magnetic susceptibility \ ~ T l l 2 . Due to the intervention of superconduc- 
tivity, it is hard to test these power laws accurately. One can deduce the continuation to 
the T-dependence below T c of C V (T) by invoking conservation of entropy on the available 
measurements and the T = limit. C V (T) ~ T 3 / 2 for T << T x is then not inconsistent 
while C V (T) ~ T clearly is inconsisitent. The temperature dependence of C v in the pseudo- 
gap regime implies that the temperature dependence and magnitude of superfluid density 
p s in the pseudogap regime are quite different from usual. This is a subject for further 
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FIG. 7: Calculated specific heat coefficient 7(T) for various dopings x compared to the values 
deduced from experiments; The specific heat data was extracted from Fig 12 of J.W. Loram, J. 
Luo, J.R. Cooper, W.Y. Liang and J.L. Tallon, J. Phys. Chem. Solids 62, 59 (2001).)- 



investigations. 



VII. CONCLUDING REMARKS 



This work rests on the development of two novel ideas using very elementary methods: 
Time-reversal violating states in a metal due to large finite range interactions and that such 
states necessarily have anisotropic gaps at the chemical potential. It is suggested that these 
determine the properties of region (II) of the phase diagram ending at the putative QCP, 
which is the basis for the singular fluctuations responsible for the properties in region (I) 
and the glue for superconductive pairing. The claim that this investigation provides the 
basic framework for a microscopic theory of the cuprates can only be made if the unusual 
predictions made are further confirmed and verified. I also discuss in this section, the 
shortcomings of the mean-field theory in this paper in relation to the properties of the 
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FIG. 8: Calculated magnetic susceptibility as a function of temperature compared with that de- 
duced from Knight shift measurements reported in Fig 3 of G. V. M. Williams, J. L. Tallon, J. 
W. Quilty, H. J. Trodahl, and N. E. Flower, Phys. Rev. Lett. 80, 377 (1998) . Note that Knight 
shift measures both the diamagnetic and the paramagnetic susceptibility. In the comparison with 
the experiments, a temperature independent diamagnetic term has been added to the calculated 
values. 

Cuprates. 

A. Predictions 

Two classes of verifiable predictions have been made in this and related work: a) The 
pseudogap state violates Time-reversal without altering the translational symmetry and ends 
at a QCP inside the superconducting dome of compositions, and b) The specific features of 
the dispersion and the resulting thermodynamic properties in the pseduogap state. 
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Concerning a), the most dramatic prediction, as has been mentioned, has been verified 28 
through specially designed ARPES experiments with polarized photons. The symmetry of 
the TRV state is consistent with the 0/j state. An important part of the verified prediction 
is that the magnitude of the measured time-reversal violation in a given direction in the 
Brillouin zone is the same for all states in the conduction band measured 2 ^ (in an energy 
range of about 0.4 eV below the chemical potential). According to the theory, these should 
all be given by the parameter 9. 

These results need to be confirmed independently with other technique, preferably more 
direct; also more compounds with a range of compositions need to be measured. Exper- 
iments using other technique a 29 ^ 30 have been proposed to test for time-reversal invariance 
and the specific point-group symmetry proposed for the current patterns. An especially 
notable prediction^ , is the nature of the superconducting state for compositions in which 
the pseudogap occurs in the normal state. This can be verified through experiments related 
to the Josephson effect. 

Several confirmations of the prediction that the pseudogap state ends at a QCP inside 
the superconducting dome now exist through tunneling 9 - and transpor t 10156 experiments in 
a magnetic field. 

Concerning b), the principal predictions are (i) that the dispersion of quasi-particles in 
the pseudogap state near the gap-edge has the form given by Eq.(70) and exhibited in 
fig. (3). To check this the spectral function including the self-energy must to be used to 
deduce the quasi-particle dispersion, (ii) that the dispersion has point-group symmetry 
in a single -domain given by Eqs. (86,87 ) of the Appendix, (iii) that the specific heat 
and the magnetic susceptibility at very low temperatures in the pseduogap state have the 
temperature dependence discussed above. The specific heat ( and the related experiment 
of thermal conductivity ) need to be done in a high enough magnetic field to suppress 
superconductivity. 

A very important issue is the specific heat at the onset of the pseudogap. As discussed 
below the transition may not have a singularity in the specific heat. But there is no reason 
why there should not be a bump. The deduction of the electronic specific heat through 
subtraction of the lattice heat capacity by comparison between samples for temperatures 
above 40-fT (where it is less than 1% of the total ) is not accurate enough to decipher bumps. 
Measurements at much lower temperatures are required which must be done in a large 
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enough field to suppress superconductivity with the transition to the pseudogap preserved. 

It is obvious that the fluctuation spectra of the form of Eq.(l) observed in the region 
(I) of the phase diagram at long wavelengths by Raman scattering should be modified in 
the classical fluctuation regime qualitatively as in Eq.(29). This is consistent with some 
recent observations^. The detailed behavior can only be predicted by a proper theory in 
the fluctuation regime (see below) and not be a mean-field theory supplemented by Gaussian 
fluctuations as in this paper. The result derived here that the fluctuations occur in specific 
IR's given by the coupling of the fluctuation vertex to the Fermions is likely to survive since 
it is based on symmetry only. 

B. Specific heat near the Transition 

A difficulty in regarding the pseduogap region as a distinct phase is that there is no 
singularity in the specific heat at the temperature where it begins to be observed in many 
other experiments. In this connection, it is worthwhile looking at the statistical mechanics for 
a TRV transition into the 0/j phase. This phase breaks both time-reversal and inversion. A 
statistical mechanical model describing the transition to such a phase is that of two-coupled 
Ising model. A realization of the model is the (anisotropic) Ashkin- Teller model: 

H = - J [(ViVj + x ^j) + yiWj^j)} (78) 

<ij> 

where the sum covers nearest neighbors on a lattice. For — 1 < y < 1 at x — 1, the model in 
two- dimensions has a line of critical points with variable critical exponents. We may identify 
(/i) as the time-reversal order parameter and (A) as the inversion breaking order parameter. 
J,x,y are positive in our problem signifying ferromangetic alignment of (/i), (A) as well as 
of (/iA). In this case the specific heat exponent a varies from to —2/3 as y varies from 
to 1, i.e. the specific heat has no divergence^. The anisotropic model (x ^ 1) also shares 
some of these special features. 

The specific heat is part of the general problem of deriving the fluctuation spectra at the 
TRV/pseudogap transition in the classical regime as well as in the quasi-classical regime, 
(i.e. Region I in the phase diagram of fig.(l)), which lead to the MFL properties as well as 
provide the glue for pairing. In this connection, it is worth noting that the Ashkin- Teller 
model has interesting topological excitations. These may well be relevant to deriving Eq. 
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(1). 



C. Shortcomings 

Some major blank spaces in the picture have already been mentioned; some others ought 
to be mentioned. The most important of these is the derivation of the phenomenolgical 
Eq.(l) and superconductivity. The fluctuation spectra of Eq. (1) has the right energy scale 
and coupling constant (both deduced from expriments) to give the right scale of supercon- 
ducting transition temperatures. The developments in this paper identify the nature of 
the fluctuations to be time-reversal and chirality fluctuations which condense to give the 
'pseudogap' region of the phase diagram. 

The major issue of principle about superconductivity is how a (nearly) momentum in- 
dependent fluctuation spectra, which predicts the observed momentum independent single- 
particle self-energy, can lead to d-wave superconductivity. The answer lies in the momentum 
dependence of the coupling functions of the current fluctuations to the fermions. As noted 
earlier, the spectra of Eq.(l) acquires a dependence on the internal momentum k of the 
particle-holes to represent the IR of the fluctuation spectra. In a separate paper, I show 
that this leads to a momentum independent self-energy but a momentum dependent pairing 
vertex. 

Some other issues worth investigating are the competition between the pseudogap state 
investigated here and the AFM state near half-filling, the transport properties in the pseu- 
dogap state and the derivation of the unusual superfluid density in the pseudogap state. A 
complete solution of the self-consistency equations in section IV-B is also desirable. Many 
of these issues which involve detailed calculations are worth the effort only if the principal 
results of this paper are verified by experiments. 
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for Condensed matter Physics at University of Karlsruhe as a Humboldt prize recipient. 
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VIII. APPENDICES 



A. a: Calculation of energies 

To determine the mean-field value of the order parameter, the dependence of the energy of 
the states of the conduction band on the invariant flux $(k) = (# x (k) —9 y (k)) is calculated. 
To do this in general , the 3x3 matrix , (17) or (21) needs to be diagonalized. This leads 
to a rather messy expression. Answers with the correct symmetry may be obtained by 
calculating the energy perturbatively in t pp /t p d- Actually, this parameter may be no smaller 
than about 1/2. Treating it as small has the advantage of obtaining results analytically and 
no disadvantage as far as the qualitative validity of the results is concerned. The change in 
energy of the conduction band is 

— 2 / 2 2 \ 2 

Se ck « 8 %^J^ cos(2(^(k) - 9 y (k)). (79) 

^ck ^xy 

Here e° k = 2i p dS xy (k) is the difference in energy of the conduction band and the non-bonding 
band for t pp = 0. x ^ y (k) for the 6/ and 6// phases are given by Eqs. (18) and (19). Riji 
and 4>i,ii are to be determined variationally. 

We write here the mean-field energy which is dependent on the order parameter up to 
order (j^) 4 - For the 6/ phase 

8F{ Rl Ai) « f - 16 £(1 - /(e ck )) S^) 2 (sin 2 7 ) _ ( 80 ) 

V k 6 ck Zt pd 




For the 9// phase, 

5F{R n An) « ^ - 16 £(1 - /(e ck )) %[&) 2 (sin 2 <f> n ) - (81) 
V k e ck ztp d 

/ Rll \i( s x c y i C x S y ) / . 2 / oj. \l SxS y ( _l \ 2 

2 (^5tH 2 ^( sm <j>IlCOs2(f) II )}^^(s x C y ±C x Sy) . 

Ztpd S xy S xy 

Compare the coefficients of the (R/t p d) 2 terms in Eqs. (80) and (81). Only the last factors 
in the parenthesis are different and the rest of the term is negative at each point (k x ,k y ). 
Then it is easy to see that the decrease in energy of the state 6// is always lower than of the 
6/ phase such that for fixed parameters the transition temperature to the former is larger 
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than the latter. Although this is true in the simplest model interactions, it is may not be 
generally true. 

On minimizing, O is ±7r/2 for either phase . Also R ^ at T = for the Q It phase 
only if: 

|< E ^#7f^(^±^) 2 (82) 

V k>k y e k S xy\ K ) 

For e d 7^ 0, tl d s 2 xy is replaced approximately by t^s^ + e 2 d . 

For 0/ phase the same condition gives a minimum value of V larger by a numerical factor 
of 0(1) with the other parameters fixed. 

Note that the x dependence occurs through that of t pd and i pp besides the Fermi-energy. 
Eq.(82) is satisfied only for x, the deviation of doping from hal-filling, less than a critical 
doping \x c . This may be estimated from the integral in Eq.(82); the critical doping at T = 
is 

x c ^VPj(2P pd ±VP pp ). (83) 

Here the plus sign is for the deviation of electron density and minus sign for the deviation 
of hole density from half-filling. 

Similarly the value of the order parameter at T = may be estimated from the energy 
expanded to order (R^) to be 

9 =< R> /2i pd ^ x c \x c - x\ 1/2 (84) 

The transition temperature T c for a fixed x can be estimated from Eqs.(80,81). The 
condition obtained can be written in terms of \x c — x\ and is given in Eq. (24). 

The numerical coeflicents in all the equations above are quite approximately determined. 
The aim here is to show only their dependence on the parameters of the model. 

B. Energy of Fluctuations 

The energy of uniform fluctuations of <fi about O — ±tt/2 can be calculated from 
Eq. (80,81) and is given, for example in the 6/ phase by 

k b xy 
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which is of O(^). 

In a given domain of the Ojj phase, with reflection symmetry broken about either k x = 
or k y = planes, the characteristic energy of oscillations to the other domain may also be 
estimated similarly. In Sec. (IV), I have called such fluctuations the mixing fluctuations. 
Their energy at long wavelengths can be calculated by allowing the order parameter to be 
a mixture of the two domains in the mean-field calculations and picking up the coefficient 
of the quadratic term (58 m ) 2 in the mixing: 

sin$(k) = %(c x /s x + c y /s y ) + 59 m (c x /s x - c y /s y ). (86) 

The fluctuation energy proportional to (S6 m ) 2 is 

Vo,m « ?A 2 E 8 (4)~ lf l^ (sin 2 ((k x - k y )a/2) cos 2 ((k x + k y )a)) (87) 

k S xy 

This fluctuation energy is seen to be smaller smaller than that of the mode of the phase of 
Time-reversal, 5<p by about a factor of 2. The fluctuation energy for the amplitude mode 
can be shown to be the largest of the frequencies of the three modes. 
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